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The unified dark energy and dark matter model within the framework of a model of a continuous 
medium with bulk viscosity (dark fluid) is considered. It is supposed that the bulk viscosity coef- 
ficient is an arbitrary function of the Hubble parameter. The choice of this function is carried out 
under the requirement to satisfy the observational data from recombination {z ~ 1000) till present 
time. 
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I. INTRODUCTION 



Modelling of an accelerated expansion of the present Universe lies on the way of creation of phenomenological 
models which may explain the observational data on one set of parameters and compare them with predictions of 
_ the models on other set. For example, a theoretical model adapts for the correct description of the acceleration of 
II i the Universe in an accessible interval z. Further, the results of modelling are extrapolated for large z which are not 
O ' accessible for observations yet. The corresponding cosmological scenario defines growth of the large-scale structure 
^ which determines the present-day fluctuations of the microwave background radiation. Certainly, the models under 
^ consideration should not contradict the available observational data within the framework of general relativity in a 
I \ field of its applicability. For the specified purposes a number of cosmological models successfully applied in the past in 
the theory of the early Universe is used. These are cosmological models with various scalar and non-scalar fields filling 
the space together with cold dark matter (see, e.g., the reviews [1]). Cosmological models of thepresent accelerated 
^ , Universe within the framework of high-order theories of gravity (HOTG) are also quite popular 

A number of models have recently been suggested [1, 0, Q which describe the present Universe with use of models 
' of a continuous medium in the presence of bulk viscosity. Consideration of effects of viscosity within the framework of 
HOTG was also carried out Note that such models were well-known in the theory of the early Universe (see, for 
[ example, 0, [I])- In particular, in Ref. Q a few exact solutions with the constant bulk viscosrh^ coefficient and with 
, the bulk viscosity being an arbitrary power function of energy density were obtained. In Ref. [3| the model of viscous 
J — ■ dark fluid is considered. The main result of this paper is the model with the constant bulk viscosity coefficient. The 
model fits the observational data on luminosity at an acceptable level. In Ref. 4] the models both with the constant 
' bulk viscosity coefficient and the bulk viscosity linearly proportional to the Hubble parameter are examined. The 
question about infiuence of viscosity on presence of a singularity in the model in the future (the so-called Big Rip) is 
investigated. 

5— i In this paper we consider a model of "viscous dark fluid" with the bulk viscosity coefficient ^.{H) which depends 

on the Hubble parameter arbitrarily. Unlike Ref. ^1, comparison of the model with the observational data is not 
restricted to the observational data on luminosity. The model is being compared with results of observations on 
change of the deceleration parameter q and values of the Hubble parameter in the range 2 > z > Q. It will be shown 
below that the model with the constant bulk viscosity coefficient does not provide a good description for q{z) and 
H{z) which follow from the observations. We propose such a dependence fi{H) which is adequate to the mentioned 
observations. The proposed model is extrapolated for z beyond the specified range 2 > z > 0. 
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II. EQUATIONS AND SOLUTIONS 



The metric of the flat Universe is taken as: 

ds^ = c^dt^ - a{tf(dx^ + dy^ + dz^). (1) 
The corresponding 0-0 component of the Einstein equations is: 

1 fdaV 8ttG 



where H is the Hubble parameter, e is the energy density of matter. Introducing the dimensionless energy density 

£* 

where the critical density = 3c^Hq /SttG (the subscript indicates the value of the parameter in the present time), 
one has from 

5=^ = h\ (3) 

Here the dimensionless Hubble parameter h is expressed in units of its present value iJp- 

The corresponding energy conservation law for the viscous dark fluid can be obtained from the equation: 

[T;' + t^1;, = 0, (4) 

where the energy-momentum tensor of matter is 
and the tensor of viscosity is 

with the bulk viscosity coefficient /i. By carrying out covariant differentiation in ([4]) with taking into account the 
metric ([T]), one can obtain the following equation: 

^+3M = 9A/i2. (5) 
do 

Here the dimensionless time 9 = Hat is introduced, and the bulk viscosity coefficient redefined with help of the 
dimensionless parameter A as 

fi = s^X/Hq. 

Similarly to let us suppose that the viscous medium has the pressure p = 0. Then equations ^ and (O imply 
one equation for the dimensionless Hubble parameter h. The case A = const is equivalent to the Murphy's model 
0. As a matter of fact this model was used in [3| for comparison with the observations. The models of the present 
Universe mentioned in Introduction use the fact that the deceleration parameter q in the past at z 3> 1 was close 
to 5 w 0.5 (cold dark matter), and in the present time the Universe expands with acceleration. For this reason all 
the models (both with scalar fields and in HOTG) are created in such a way that the present inflation appears for 
rather small values of the Hubble parameter, or what is the same thing, for small average density of matter. In 
this connection we will consider the model ([5]) in which the dimensionless bulk viscosity of the dark fluid A is not a 
constant but an arbitrary function of the parameter h. In this paper we chose this function as: 

9A = 3tanhf-^Y (6) 



where 6, n are arbitrary constants which will be defined from the observational data later. 
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For comparison with the observational data, it is convenient to rewrite equations ([3]) and ([5|) through the redshift 
z — 1/a ^ 1. Then one has the equation for h in the form: 



2(z + 1)^ + 3/i = 3tanh 
dz 



and the deceleration parameter q will be: 



{z + 1) dh 
h dz 



(7) 



The parameters b, n in the model are being chosen from the requirement that in the present time (at h = 1 with 
account of selected normalization) the deceleration parameter q should be close to the observable value q « —0.6. 
This value of q can be obtained at & = 0.95 and n — 2. Using these parameters, we have compared the model 
under consideration with the observational data in the range 2 > z > 0. The results are presented in Fig. For 
comparison with our model, the ACDM model with the same "initial conditions" is also shown in Fig. ([T]). The model 
with the constant bulk viscosity from [sj is presented as well. 
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FIG. 1: The deceleration parameter for our model (thick solid line), the model with the constant bulk viscosity from 0] 
(dashed line), the ACDM model (dashed-dot line). The central solid thin line represents the best-fit, the light grey contours 
represent the Icr confidence level, and the dark grey contours represent the 2a confidence levels (the data taken from ^0|). 

From Fig. ([T]), it follows that the model with the constant bulk viscosity deviates from the observations appreciably. 
The model under consideration with \(K) lies within the confidence levels. Note that this model close to the ACDM 
model (with fi^ = 0.3 and JIa = 0.7), although there is no any special A-term in the model. Outside the interval of 
observations (in practice at z > 2) our model, the ACDM model and extrapolation of the observable data from Q 
are close to each other and give q w 0.5. The model with constant A has much less value of q and reaches q « 0.5 at 
z > 40. 

An extrapolation in the future for all three models is also shown in Fig. ([T]). Asymptotically (z ^ — 1) all three 
models tend to the de Sitter model. 



III. CONCLUSION 



The bulk viscosity in our model is an example of a dynamic A-term. However, our model is close to the ACDM 
model. As it was rightly noted in [sf, the model with the viscosity does not give possibility to divide the true dust 
filling the Universe, and dark matter generated by the bulk viscosity. That is why it is difficult to introduce a 
phenomenological equation of state p — we which is often used for interpretation of the observational data. Influence 
of the bulk viscosity on formation of the large-scale structure of the Universe demands special examination. But 
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taking into account that the viscosity is being "involved" at rather small z « 2, it is possible to expect its influence 
on the dynamics of galactic clusters only at later non-linear stage. 
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